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ABSTRACT 

The intense, bursting X-ray source in the globular cluster NGC 1851 was one of the first cluster sources 
discovered, but has remained optically unidentified for 25 years. We report here on results from Hubble Space 
Telescope WFPC2 multicolor images in NGC 1851. Our high spatial resolution images resolve ~200 objects in the 
3" radius Einstein X-ray error circle, 40 times as many as in previous ground-based work. A color-magnitude 

diagram of the cluster clearly reveals a markedly UV-excess object with B ~ 21, (U - B) 0.9, only 2" from 

the X-ray position. 

The UV-excess candidate is 07 1 2 distant from a second, unremarkable star that is 0.5 mag brighter in B\ thus 
ground-based studies of this field are probably impractical. Three other UV-excess objects are also present 
among the "-16,000 objects in the surveyed region of the cluster, leaving a ~5% probability that a UV-excess 
object has fallen in the X-ray error circle by chance. No variability of the candidate is seen in these data, although 
a more complete study is required. If this object is in fact the counterpart of the X-ray source, previous inferences 
that some globular cluster X-ray sources are optically subluminous with respect to low-mass X-ray binaries in the 
field are now strengthened. 

Subject headings: globular clusters: individual (NGC 1851) — stars: neutron — ultraviolet: stars — 

X-rays: bursts — X-rays: stars 


1. INTRODUCTION 

It has been known for over two decades that cluster X-ray 
sources are overabundant with respect to those in the field 
(Katz 1975; Clark 1975). While globular clusters only contain 
10 4 of the mass of the Galaxy, over 10 2 * of the bright X-ray 
sources (L x > 10 26 ergs s 1 ) are found in globulars. The cause 
of this discrepancy may be a difference in the nature of 
globular cluster X-ray sources from those in the field or more 
likely an enhancement of their formation probability. 

Despite a quarter-century of effort, there are very few 
intense globular cluster X-ray sources with convincingly iden- 
tified optical counterparts. AC 21 1 in M15 (Auriere, Le Fevre, 
& Terzan 1984; Charles, Jones, & Naylor 1986) was discovered 
over a decade ago; it is surprisingly optically luminous. A 
variety of candidates in NGC 6712 were offered (e.g., Cud- 
worth 1988; Nieto et al. 1990; Bailyn et al. 1991; Auriere & 
Koch-Miramond 1992) and the counterpart, star S, was even- 
tually unambiguously confirmed by Anderson et al. ( 1 993) with 
HST images. A counterpart for the X-ray source in NGC 
6624 was discovered to dominate the cluster light at 1400 A in 
HST Faint Object Camera images by King et al. (1993). These 
three optical counterparts discovered so far show a wide range 
in optical luminosity, and recently the spectra of AC 21 1 and 
star S were found to differ significantly (Downes et al. 1996). 

The scarcity of optical counterparts is almost certainly 

1 Based on observations with the NASA/ESA Hubble Space Telescope 

obtained at the Space Telescope Science Institute, which is operated by the 

Association of Universities for Research in Astronomy, Inc., under NASA 

contract NAS5-26555. 


primarily due to the extreme crowding in the clusters, which 
limits the utility of ground-based identification programs. As 
part of the Hubble Space Telescope Faint Object Spectrograph 
team’s observing program, we have obtained multicolor HST 
images of several globular clusters containing X-ray sources, 
with the aim of resolving optical counterparts in crowded fields 
and studying their spectra. We report here on the results from 
WFPC2 images in NGC 1851. 

X-ray observations of the strong source in NGC 1851 are 
virtually as old as satellite X-ray astronomy, dating back to the 
days of OSO 7 (Clark, Markert, & Li 1975); Forman & Jones 
(1976) first observed the bursting nature of X0512-401 with 
Uhuni. Recent X-ray observations of the system have been 
made and reviewed by Callanan, Penny, & Charles (1995). 
This source would seem a particularly good candidate for an 
optical identification, as there is an accurate Einstein HRI 
X-ray position (Grindlay et al. 1984), and it lies over 2 core 
radii (M2") north of the cluster center. 

Ground-based imaging studies in NGC 1851 by Auriere, 
Bonnet-Bidaud, & Koch-Miramond (1994, hereafter ABK) do 
not reveal any especially good optical candidate (e.g., no 
object with marked UV excess) for the X-ray source 
X0512-401. However, these authors call attention to an 
object they denote as XI, as the best candidate they are able 
to resolve in a 3" error circle. As its color (U - B 0) is no 
more UV than many hot blue horizontal-branch stars in the 
cluster, ABK speculate that XI may be a new sort of low-mass 
X-ray binary (LMXRB), in which there is very little X-ray 
heating of the accretion disk. 

Previous studies of NGC 1851 using pre-servicing-mission 
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HST images did not detect any better candidates (Margon et 
ai. 1992; Deutsch et al. 1994). The WF/PC images used in 
those studies do not go nearly as deep as the WFPC2 data used 
in the present work. Deutsch et al. (1994) reported several 
objects that were bluer than XI at the detection threshold of 
those images, but were unable to isolate any single object 
unusual enough to be considered a better candidate. These 
earlier studies would have easily revealed counterparts as 
luminous as AC 211 or star S. Clearly, deeper exposures are 
required to look for fainter objects in NGC 1851. 

2. OBSERVATIONS AND DATA REDUCTION 

2. 1 . Planetary Camera Imagery 

On 1996 April 10, we obtained HST WFPC2 images of 
NGC 1851. Two 900 s exposures through the F336W filter 
were taken, followed by four 300 s F439W exposures, and 
finally two more 900 s F336W exposures; the F336W and 
F439W filters are similar to Johnson U and B , respectively. 
The frames have been processed through the standard data 
reduction pipeline at STScI. Further reduction was performed 
with software written in interactive data language (IDE) by 
E. W. D. or available in the IDL Astronomy User’s Library 
(Landsman 1993). Some photometry was performed with 
DoPHOT (Schechter, Mateo, & Saha 1993). 

Although the images were obtained with the fine-lock mode 
of HST, the point-spread functions (PSFs) clearly show elon- 
gation in the Y direction in the Planetary Camera (PC) images, 
with major- to minor-axis ratios of about 1.25. Furthermore, 
while the images show excellent registration in theX direction, 
Y offsets between exposures through the same filter are as high 
as 0.4 pixels. 

For each filter, the set of four images is combined into one 
master image using software written by E. W. D. First, 
cosmic-ray hits are identified in each image, and those pixels 
are set to zero in the data and “exposure mask” arrays. The 
identification algorithm is similar to the one described in Saha 
et al. (1996), and is carefully adjusted so that the under- 
sampled cores of bright stars with subpixel shifts between 
images are not mistaken for cosmic-ray events. The four 
images and masks are aligned to ~0.01 pixel accuracy and 
summed. The summed frame is then divided by the mask to 
obtain the final image. In this way, cosmic-ray event pixels are 
excluded so that the information in neighboring pixels is not 
contaminated during interpolation. 

2.2. Astrometry 

In order to determine the position of the X-ray source 
coordinates on the HST PC image, we establish a coordinate 
system based on the HST Guide Star Catalog (GSC) reference 
frame. We begin with the astrometric solution from a digitized 
Quick V image used to generate the HST Guide Star Catalog 
(Lasker et al. 1990); this image, obtained directly from STScI, 
contains the astrometric solution in the image header. By 
centroiding 29 stars in the Quick V image and the correspond- 
ing objects in a short ground-based /?-band CCD image, kindly 
provided by H. Ford, the astrometric reference frame is 
transferred to the CCD image with an error in the solution of 
OH. Next, the astrometric frame is transferred from the CCD 
image to the PC image using a set of nine common isolated 
stars, with error in the solution of 0"02, negligible compared to 
the first step. The final result is that coordinates may be 
determined on the PC image to within O'.'I in the GSC 


reference frame. However, there may well be some systematic 
offset, a ~ 0"5, from frames based on other reference catalogs 
(Russell et al. 1990). 

We calculate the weighted mean of the X-ray coordinates 
from reprocessed Einstein HRI data in the HRICFA database 
obtained through the High Energy Astrophysics Science Ar- 
chive Research Center Online Service, provided by the NASA 
Goddard Space Flight Center. These coordinates are not 
significantly different from those in Grindlay et al. (1984), so 
we use the published coordinates in this work. The 90% 
confidence radius is reported by Grindlay et al. (1984) as 3". 

The optical position of the X-ray coordinates chosen by 
ABK, which we estimate from Figure 1 of that paper, is ~1" 
northeast of our location, which is not unreasonable given the 
different source of astrometric reference. 

The left panel of Figure 1 (Plate Lll) shows a 24" X 24" 
region of the cluster through the F336W filter. A 3" radius 
error circle is overlaid, using our astrometric solution. The 
right panel of Figure 1 shows a 8" X 8" region of the cluster 
through the F336W filter, centered on our optical position for 
the X-ray source, which we have marked with a cross. 

2.3. Photometry 

Photometry was performed on the combined F336W and 
F439W PC frames with the DoPHOT software (Schechter et 
al. 1993) and a PSF function-fitting procedure written by 
E. W. D. The results are similar, although DoPHOT was not 
able to converge accurately on some of the brighter objects; 
however, it performed better on the fainter objects. For this 
reason, we use the results derived from DoPHOT. The fitted 
magnitudes were calibrated with aperture photometry of sev- 
eral isolated objects. Aperture corrections are taken from 
Table 2a in Holtzman et al. (1995b). The photometric mea- 
surements have not been corrected for geometric distortions in 
the PC, but the simple correction for charge transfer efficiency 
losses detailed in Holtzman et al. (1995b) has been applied. 
We use the photometric zero points for the STMAG system 
from Table 9 in Holtzman et al. (1995a). Systematic errors for 
all magnitudes due to uncertainties in detector performance, 
absolute calibration, and filter transformations are ~5%. 

The photometric results are displayed in Figure 2 as a 
color-magnitude diagram (CMD) of all ~5500 objects with 
reliable measurements in the PC frame. Objects which are 
mentioned in the discussion below are marked in the diagram. 

3. DISCUSSION 

While ABK were only able to resolve five stars within a 3" 
radius error circle, we are able to resolve and measure 
magnitudes for 194 objects with 15 ^ ^439 < 23. 

Star XI, suggested by ABK as the best candidate, does have 
the brightest m 336 magnitude within the error circle. However, 
its color and magnitude are unremarkable compared with 
other cluster horizontal-branch stars. We resolve the star 
labeled as X2 by ABK into two main components, which we 
label X2a and X2b. Additional, fainter stars also contribute 
light to X2. Object X2a is nearly in the center of our error 
circle, but has normal colors. In the CMD, X2b is found to 
have a more UV color than star XI, but is also unremarkable 
compared with other stars in the cluster (Fig. 2). 

In the CMD of the entire 34" X 34" PC field, there are two 
objects that stand out as having marked UV excess, with colors 
similar to known LMXRBs. One of these is 2715 away from 
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Fit;. 2. — Color-magnitude diagram for ~5500 stars in the PC frame. 
Magnitudes are in the STMAG system; see text for details. All objects which 
are discussed in the text are labeled here, including star A, which we select as 
the optical counterpart due to a significant UV excess. An approximate \f [{ 
scale is also provided, assuming (m M ) ~ 15.43, E(B — V) = 0.02, and 
(B - mw) = 0.65. This latter correction changes slightly with stellar color; 
0.65 is appropriate for F type through the hottest stars, while 0.4 is more 
appropriate for MO stars. We also add the distance- and reddening-corrected 
magnitude and color of star S in NGC 6712 for comparison. 

our optical position for the X-ray coordinates; this is well 
within the 90% confidence error circle of the X-ray coordi- 
nates. We denote this object as star A and provide magnitude 
measurements in Table 1 . Figure 3 (Plate LI 1) shows a 3" X 3" 
region of our PC data centered on star A. We also label star B, 
which is probably unrelated to, but blended with, star A, at a 
separation of only O' T2. The photometry of these objects is 
slightly hampered by the fact that the elongation of PSFs, 
discussed in §2.1, is in nearly the same direction as the 
orientation of these two stars. Future observations will be 
complicated by the proximity of the two objects. 

We have also examined the data for possible time variability 
of Star A. The first and last F336W exposures are separated by 
3 hr, while the first and last F439W exposure are separated by 
only 70 minutes. Neither series shows any significant evidence 
for variability, although these data do not allow an adequate 
study. The 1 a photometric uncertainty from photon counting 
statistics is only 2% and 3% for star A for the F336W and 
F439W filter exposures, respectively. However, since star A is 
blended with star B, the measurements have a higher uncer- 
tainty. On the final F336W exposure, star A was contaminated 
with a cosmic-ray event, increasing the error of the measure- 
ment for that exposure. We set a variability upper limit of 
0.10 mag during the time of our observations. These data 
could not have detected any small-amplitude, short-term 
variability of the scale reported by Homer et al. (1996) for the 
globular cluster X-ray source counterpart star S in NGC 6712. 
Archival WF/PC images of NGC 1851 have also been exam- 
ined for variability. Star A is above the detection threshold 
only in F336W exposures taken on 1993 August 21 . There is no 
obvious evidence that star A has changed in brightness during 
this 2.6 yr interval, but the variability upper limit is not very 
stringent, ~0.4 mag. 


TABLE 1 

Photometry and Coordinates for Selected Obje cts in NGC 1851 

ct( 2t)(X)) 8(2000) 

5 h 14 m -40°02 f 


Object (s) (arcsec) ~ w 4 .w) 


Star A a 6.428 38.38 10.73 20.46 -0.73 

StarB 6.425 38.26 20.43 19.92 0.51 

Star Xl b 6.523 35.81 16.58 15.52 1.06 

Star X2a c 6.556 36.81 17.51 16.40 1.11 

Star X2b c 6.592 36.59 16.85 15.92 0.93 


Optical counterpart for X05 12-401 suggested in this work. 

b Designation from Auriere el al. 1994. 

1 X2, a designation from Auriere et al. 1994, is multiply resolved in the 
current work. 

Since star A is not the only UV-excess object in its part of 
the CMD, there is some probability that it is not the X-ray 
counterpart but rather an unrelated object similar to the other 
cluster UV objects. To estimate the probability, we measure 
magnitudes for all stars in the PC field as well as the three WF 
fields. Although this larger sample has 13 times more area, it 
contains only 3 times more stars than the PC field alone, as the 
latter is centered near the cluster core and there is a sharp 
radial gradient of stellar density. Of the 16,000 stars for which 
we obtain good measurements, we find four objects with 
significant UV excess of “ w 4 v>) ^ —0.5 in the range 
15 < m 4 M ^ 22. We also search each of the fields with an 
image-subtraction technique to further ensure that no UV- 
excess objects are missed. 

Therefore, approximately 2.5 x 10 4 stars with m 43 y < 22 in 
the cluster are markedly UV excess, and as we measure 
194 stars within the error circle, we estimate a ~5% probabil- 
ity that one of these UV-excess objects unrelated to the X-ray 
source will be found in the error circle by chance. The hazards 
of these a posteriori probability estimates are well known, so 
this result must surely be regarded as qualitative rather than 
quantitative. The nature of the faint UV-excess objects ob- 
served by //.ST in other globular cluster cores has been recently 
discussed by numerous authors. 

In NGC 6712 the markedly UV-excess optical counterpart 
star S is surprisingly optically subluminous as compared with 
LMXRBs in the field (Anderson et al. 1993). If the magnitude 
and color of star S reported in Homer et al. (1996) are 
corrected for the difference in reddening and distance 2 be- 
tween NGC 6712 and NGC 1851, one estimates that star S, 
were it present in NGC 1851, would have w 4 .™ = 19.4, 
(mu,, - m 4 ») = -0.9 in this CMD. Therefore, while the 
colors of star S and the candidate presented here are similar, 
star S is more luminous by 1 mag. If star A is the correct 
identification of the NGC 1851 X-ray source, the luminosity 
discrepancy of this subclass of objects from the field LMXRBs 
is further exacerbated. 

To quantify this issue, we use the STSDAS SYNPHOT 
package and find that the blackbody spectrum which best fits 
the star A color {nu M) - w 4>) ) = -0.73 has T ~ 20,000 K. 
Using a thermal spectrum of this temperature, we derive 
approximate relations to convert to Johnson U and B magni- 
tudes: (U - m.w,) = 0.43 and (B - m w ) = 0.65. Applying 
these zero-point offsets, we estimate V = 20.2, B = 21.1, and 
(U — B) = -0.9. A typical LMXRB in the field has 

- A distance of 12.2 kpc and E(B - 1 ) = 0.02 for NGC 1851, and 6.8 kpc 
and E(B V) = 0.46 for NGC' 6712. were adopted from Peterson { 1993). The 
reddening curve of Savage & Mathis (1979) was used. 
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My ~ 1 .2, (B - V) ~ 0, and (£/ - B) ~ -0.9 (van Paradijs 
1983). Therefore, the NGC 1851 candidate has ( U - B) 
typical of LMXRBs but is underluminous by over 4 mag, with 
M b = 5.6. One model for the system described by Callanan et 
al. (1995) invokes a subgiant secondary star; our faint coun- 
terpart rules out this possibility. 

4. CONCLUSION 

We recompute and confirm the previously published optical 
position of the Einstein HRI 3" radius error circle for the 
bright X-ray source in NGC 1851 using the HST GSC as a 
reference frame. With WFPC2 observations, 194 stars are 
resolved and photometered within the error circle. Star XI, 
the best candidate previously offered in the literature, is 
resolved as a single star, but is not remarkable in a CMD of the 
region. Star X2 is resolved as two separate stars, the brighter 
of which is more UV-excess than XI but still is not remark- 
able. 

However, we find two UV-excess objects in the PC frame, 
one of which is only 2'.T5 from the X-ray source position. 
Therefore, we suggest this as a much more likely candidate for 
the optical counterpart of the X-ray source, and we denote it 
as star A. We search for variability, and find none, for this 
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PLATE LI 1 



Fig. 1 .—Left panel : 24" x 24" section of the HST PC image of NGC 1851 with the F336W filter. The image has been reoriented so that north is up and east is 
to the left. A 3" radius error circle is drawn about the position of the X-ray coordinates. Right panel: 8" X 8" section of the same image centered at the X-ray 
coordinates. A cross denotes the optical position, derived in this paper from the GSC source data, of the X-ray coordinates, and the 3" radius error circle is overlaid. 
Several objects discussed in ihe text are labeled. Star A is selected as the optical counterpart due to a significant UV excess. 

Dan sat et al. (see 472, L98) 



Fig. 3._3" x 3" sections of the F336W (U) and F439W {B) images centered on star A, which we select as the optical counterpart due to a significant UV excess. 
These images are still in the original orientation, so north and east are labeled. Star B, 07 12 distant, is apparently unremarkable, but may complicate future 
observations. 
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